CLAIMS 



1 . A method for selecting a near optimal or optimal mathematical model 
from a set of candidate models, comprising: 

defining a candidate search space having n dimensions, wherein n is a positive 
integer and each dimension represents a set of mutually exclusive features from which 
one is chosen for each candidate model. 

2. The method of claim 1, wherein said candidate search space is searched 
for a near optimal or optimal model. 

3. The method of claim 2 wherein said search is accomplished by a method 
selected from the group consisting of: full grid search, simulated annealing, integer 
programming, scatter search/path relinking, neural networks, tabu search and genetic 
algorithm. 

4. A method for automated generation of NONMEM/NMTRAN control 
files, comprising: 

a) selecting one feature from each of n sets of candidate features, wherein n is a 
positive integer; and 

b) substituting text associated with each selected feature into a control file 

template. 

5. A method for automated evaluation of the optimality of a model 
comprising: 

combining the log likelihood value with, optionally, a penalty for each parameter 
estimated, optionally, a penalty for each element of the interindividual variance matrix 
estimated, optionally, a penalty for each element of the intraindividual variance matrix 
estimated, optionally, a penalty imposed if the minimization does not conclude 
successfully, optionally, a penalty if the standard errors of the parameter estimates cannot 
be obtained, optionally, a penalty if the correlation matrix of the estimates has any 
element > 0.95, and optionally a "niche" penalty for being similar to other models in the 
population (within a "niche radius" of other models). 
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6. A method for selecting a near optimal or optimal mathematical model 
from a set of candidate models, comprising: 

a) defining a candidate search space having n dimensions, wherein n is a positive 
5 integer and each dimension represents a set of mutually exclusive features from which 

exactly is chosen for each candidate model and each model is represented by a bit string; 

b) assessing the fitness of each model in said population; 

c) optionally, scaling the fitness of each model to be between and upper limit R 
and a lower limit S wherein the ratio of R to S is between 2: 1 and 100: 1 ; 

10 d) providing a number y of models to be in a subsequent generation; 

e) selecting with replacement y number of parents of the said subsequent 
generation from the current generation, wherein the probability of selection of a model in 
the current generation is proportional to said fitness or optionally to said scaled fitness; 

f) associating said parents into m groups comprising p parents where p is an 
15 integer greater than 1 ; 

g) optionally, selecting some fraction of the m groups of parents to undergo at 
least one cross over; 

h) optionally, crossing over said selected fraction at a random location on said bit 
string to create two new individuals for said subsequent generation; 

20 i) assigning bit strings in current generation that are not selected for cross over to 

said subsequent generation; 

j) optionally, randomly mutate bits of said subsequent generation bit strings 

wherein said mutation comprises change a bit value 0 to a bit value of 1 or changing a bit 

value of 1 to a bit value of 0; and 
25 k) repeating the steps of b through j until further improvement in mean and 

maximum fitness no longer occurs. 

7. The method of claim 6, wherein said initial population is a random 
population. 



42 



Docket No: MS1US 



8. The method of claim 6, wherein said fitness is assessed by calculating 
some statistic of the goodness of fit of the model to the data and, optionally, adding cost 
associated with desirable atributes of the model, including parsimony (fewer parameters). 

5 9. The method of claim 8, wherein the goodness of fit of the model to the 

data is the log likelihood of the data, given the model. 

10. The method of claim 6 wherein the ratio of R to S is between 10: 1 and 

50:1. 

10 

11. The method of claim 6, wherein the number of models in the subsequent 
generation is equal to the number of models in the current generation. 

12. The method of claim 6 wherein p = 2. 

15 

13. The method of claim 6 wherein said fraction to undergo at least one cross 
over is selected randomly. 

14. The method of claim 6 wherein said fraction to undergo at least one cross 
20 over is between 0.4 to 1 .0 

15. The method of claim 6 wherein said models represent pharmacokinetic 
models or pharmacodynamic models. 

25 16. The method of claim 15 wherein said sets of mutually exclusive features 

comprise one or more members of the group consisting of: the number of 
pharmacokinetic compartments, the presence of non-linear elimination, the presence of 
non-linear absorption, the presence of interindividual variability on each parameter, the 
function describing the interindividual variability of each parameter, the function 

30 describing the residual variability, the structure of the interindividual covariance matrix, 
emax pharmacodynamic model, linear pharmacodynamic model, types 1 through 4 
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indirect response pharmacodynamic model, the presence of an effect compartment, the 
relationship between drug elimination and age, the relationship between drug elimination 
and renal function, the relationship between drug elimination and liver function, the 
relationship between drug elimination and gender, the relationship between drug 
5 elimination and weight, the relationship between drug volume of distribution and age, the 
relationship between drug volume of distribution and gender, the relationship between 
drug volume of distribution and weight, the relationship between drug volume of 
distribution and cardiac function, the relationship between drug volume of distribution 
and renal function, the relationship between drug volume of distribution and liver 
10 function, the relationship between drug bioavailability and age, the relationship between 
drug bioavailability and gender, the relationship between drug bioavailability and weight, 
the relationship between drug bioavailability and liver function, the relationship between 
drug bioavailability and renal function. 

15 17. A method for selecting a near optimal or optimal mathematical model 

from a set of candidate models, comprising: 

a) defining a candidate search space having n dimensions, wherein n is a positive 

integer and each dimension represents a set of mutually exclusive features from which 

exactly one is chosen for each candidate model; and 
20 b) searching the candidate search space using simulated annealing, wherein 

simulated annealing comprises the steps of: 

i) randomly selecting one model from the candidate search space; 

ii) selecting an initial value for temperature (T) wherein T represents the 
tolerance of a minimization process for retaining a change in the model 

25 that results in a higher energy; 

iii) assessing the energy of the initial model; 

iv) randomly changing the model to generate a subsequent model; 

v) assessing the energy of the subsequent model; 

vi) retaining the subsequent model as the current model if the energy is 
30 lower than the current model; 
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vii) if the energy of the subsequent model is higher than the energy of the 
current model the probability of retaining it is equal to: 

AE 

where T is the temperature, AE is the change in energy (current model 
energy - subsequent model energy), and k is Boltzmans constants, or 
5 Otherwise, rejecting the subsequent model; 

viii) reducing the value of T; and 

ix) repeating the steps of iv through viii until further reduction in energy 
no longer occurs. 

10 18. The method of claim 17 wherein said sets of mutually exclusive features 

comprise one or more members of the group consisting of: the number of 
pharmacokinetic compartments, the presence of non-linear elimination, the presence of 
non-linear absorption, the presence of interinOdividual variability on each parameter, the 
function describing the interindividual variability of each parameter, the function 

15 describing the residual variability, the structure of the interindividual covariance matrix, 
emax pharmacodynamic model, linear pharmacodynamic model, types 1 through 4 
indirect response pharmacodynamic model, the presence of an effect compartment, the 
relationship between drug elimination and age, the relationship between drug elimination 
and renal function, the relationship between drug elimination and liver function, the 

20 relationship between drug elimination and gender, the relationship between drug 

elimination and weight, the relationship between drug volume of distribution and age, the 
relationship between drug volume of distribution and gender, the relationship between 
drug volume of distribution and weight, the relationship between drug volume of 
distribution and cardiac function, the relationship between drug volume of distribution 

25 and renal function, the relationship between drug volume of distribution and liver 

function, the relationship between drug bioavailability and age, the relationship between 
drug bioavailability and gender, the relationship between drug bioavailability and weight, 
the relationship between drug bioavailability and liver function, the relationship between 
drug bioavailability and renal function. 
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19. A method for selecting a near optimal or optimal mathematical model 
from a set of candidate models, comprising: 

a) defining a candidate search space having n dimensions, wherein n is a positive 
5 integer and each dimension represents a set of mutually exclusive features from which 

exactly one is chosen for each candidate model; and 

b) searching the candidate search space using full grid search wherein full grid 
comprises the evaluation of every possible model in the search space. 

10 20. The method of claim 19 wherein said sets of mutually exclusive features 

comprise one or more members of the group consisting of: the number of 
pharmacokinetic compartments, the presence of non-linear elimination, the presence of 
non-linear absorption, the presence of interindividual variability on each parameter, the 
function describing the interindividual variability of each parameter, the function 

15 describing the residual variability, the structure of the interindividual covariance matrix, 
emax pharmacodynamic model, linear pharmacodynamic model, types 1 through 4 
indirect response pharmacodynamic model, the presence of an effect compartment, the 
relationship between drug elimination and age, the relationship between drug elimination 
and renal function, the relationship between drug elimination and liver function, the 

20 relationship between drug elimination and gender, the relationship between drug 

elimination and weight, the relationship between drug volume of distribution and age, the 
relationship between drug volume of distribution and gender, the relationship between 
drug volume of distribution and weight, the relationship between drug volume of 
distribution and cardiac function, the relationship between drug volume of distribution 

25 and renal function, the relationship between drug volume of distribution and liver 

function, the relationship between drug bioavailability and age, the relationship between 
drug bioavailability and gender, the relationship between drug bioavailability and weight, 
the relationship between drug bioavailability and liver function, the relationship between 
drug bioavailability and renal function. 

30 
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21 . A method for selecting a near optimal or optimal mathematical model 
from a set of candidate models, comprising: 

a) defining a candidate search space having n dimensions, wherein n is a positive 
integer and each dimension represents a set of mutually exclusive features from which 

5 exactly one is chosen for each candidate model; and 

b) initializing the search with a call to OCLSetup in the OptQuest callable library 
and initialize a population of models with a call to OCLInitPop; 

c) initializing each search dimension with a call to OCLDefineVar in the 
OptQuest callable library; 

10 d) selecting an initial model from the candidate search space using scatter 

search/path relinking and tabu search as implemented in the OptQuest Callable library 
from OptTek systems by calling the function OCLGetSolution; 

e) searching the candidate search space using Scatter search/path relinking/Tabu 
search using the OptQuest Callable library wherein Scatter search/path relinking/Tabu 

15 search comprises the steps of: 

i) evaluating the fitness of the current model; 

ii) adding the value of the fitness of the current model to the OptQuest 
Callable library database with a call to the function OCLPutSolution; 

iii) finding the fitness of the best model thus far evaluated with a call to 
20 the function OCLGetBest in the OptQuest Callable Library; 

iv) getting the subsequent model with a call to the function 
OCLGetSolution; and 

v) repeating steps i-iv until either the required number of evaluations or 
convergence is seen; and 

25 f) deleting current problem from memory with a call to OCLGoodBye. 

22. The method of claim 21 wherein said sets of mutually exclusive features 
comprise one or more members of the group consisting of: the number of 
pharmacokinetic compartments, the presence of non-linear elimination, the presence of 

30 non-linear absorption, the presence of interindividual variability on each parameter, the 
function describing the interindividual variability of each parameter, the function 
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describing the residual variability, the structure of the interindividual covariance matrix, 
emax pharmacodynamic model, linear pharmacodynamic model, types 1 through 4 
indirect response pharmacodynamic model, the presence of an effect compartment, the 
relationship between drug elimination and age, the relationship between drug elimination 

5 and renal function, the relationship between drug elimination and liver function, the 
relationship between drug elimination and gender, the relationship between drug 
elimination and weight, the relationship between drug volume of distribution and age, the 
relationship between drug volume of distribution and gender, the relationship between 
drug volume of distribution and weight, the relationship between drug volume of 

10 distribution and cardiac function, the relationship between drug volume of distribution 
and renal function, the relationship between drug volume of distribution and liver 
function, the relationship between drug bioavailability and age, the relationship between 
drug bioavailability and gender, the relationship between drug bioavailability and weight, 
the relationship between drug bioavailability and liver function, the relationship between 

15 drug bioavailability and renal function. 

23. A computer program product that selects a near optimal or optimal 
mathematical model from a set of candidate models, the computer program product 
comprising a computer usable storage medium having computer-readable program code 
20 embodied in the medium, the computer-readable code comprising: 

computer-readable code that is configured to accept a user defined candidate 
search space having n dimensions, wherein n is a positive integer and each dimension 
represents a set of mutually exclusive features from which exactly one is chosen for each 
candidate model. 

25 24. The method of claim 23 wherein computer-readable code is configured to 

search said space for a near optimal or optimal model 

25 The method of claim 24 wherein said search is conducted using a method 
selected from the group consisting of: full grid search, simulated annealing, integer 
programming, scatter search/path relinking, neural network, tabu search and genetic 

30 algorithm. 
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26. A computer program product that selects a near optimal or optimal 
mathematical model from a set of candidate models, the computer program product 
comprising a computer usable storage medium having computer-readable program code 
embodied in the medium, wherein the computer-readable code is configured to: 
5 a) accept user input defining a candidate search space having n dimensions, 

wherein n is a positive integer and each dimension represents a set of mutually exclusive 
features from which exactly one is chosen for each candidate model and each model is 
represented by a bit string; 

b) assess the fitness of each model in said population; 
10 c) optionally, scale the fitness of each model to be between and upper limit R and 

a lower limit S wherein the ratio of R to S is between 2: 1 and 100: 1 ; 

d) provide a number y of models to be in a subsequent generation; 

e) select with replacement y number of parents of the said subsequent generation 
from the current generation, wherein the probability of selection of a model in the current 

15 generation is proportional to said fitness or optionally to said scaled fitness; 

f) associate said parents into m groups comprising p parents where p is an integer 
greater than 1; 

g) optionally, select some fraction of the m groups of parents to undergo at least 
one cross over; 

20 h) optionally,cross over said selected fraction at a random location on said bit 

string to create two new individuals for said subsequent generation; 

i) assign bit strings in current generation that are not selected for cross over to 
said subsequent generation; 

j) optionally, randomly mutate bits of said subsequent generation bit strings 
25 wherein said mutation comprises change a bit value 0 to a bit value of 1 or changing a bit 
value of 1 to a bit value of 0; and 

k) repeat the steps of b through j until further improvement in mean and 
maximum fitness no longer occurs. 
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